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Abstract 

An abstract approximation and convergence theory for the closed-loop solution of discrete-time 
linear-quadratic regulator problems for parabolic systems with unbounded input is developed. 
Under relatively mild stabilizability and detectability assumptions, functional analytic, operator 
theoretic techniques are used to demonstrate the norm convergence of Galerkin-based 
approximations to the optimal feedback control gains. The application of the general theory to a 
class of abstract boundary control systems is considered. Two examples, one involving the 
Neumann boundary control of a one dimensonal heat equation, and the other, the vibration control 
of a cantilevered viscoelastic beam via shear input at the free end, are discussed. 
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1. Introduction 

In this paper we develop an abstract approximation framework for linear quadratic regulator 
(LQR) problems for infinite dimensional discrete-time parabolic systems with unbounded input. 
More specifically, we consider the application of the abstract approximation theory developed in 
[9] to the class of systems which are open-loop abstract parabolic and whose input operators have 
range in some space larger than the standard state space in which the problem is usually formulated. 
The theory we present here is a discrete-time analog of the results given in [3] for continuous-time 
parabolic systems. However, in contrast to the treatment in [3] which is restricted to the case of 
bounded input, we on the other hand are able to handle a relatively wide class of systems involving 
unbounded input. Our framework is applicable for example, to a variety of boundary control 
systems. 

The abstract framework we develop here leads to a set of three relatively easily verified 
conditions (stabilizability, detectability, and strong convergence of the othogonal projections 
corresponding to the approximating Galerkin subspaces) which when satisfied yield norm 
convergence of finite dimensional approximations to the optimal feedback gains. We employ a 
functional analytic, operator theoretic approach to obtain a relatively complete and reasonably 
general theory which is easily applied in practice to a wide class of important problems. 

An outline of the remainder of the paper is as follows. In section 2 we briefly outline the 
abstract theory and approximation results for infinite dimensional discrete-time LQR problems 
developed in [9]. Section 3 is concerned with abstract parabolic control systems with unbounded 
input, their operator theoretic and discrete-time formulation, and the associated infinite dimensional 
optimal control problem. The approximation and convergence theories are discussed in section 4. 
In section 5 we consider abstract boundary control systems and present two examples; one 
involves the Neumann boundary control of a one dimensional heat equation, and the second is 
concerned with the vibration control of a cantilevered viscoelastic beam via a shear input at the free 
end. In section 6 we summarize our findings and make some concluding remarks. In an appendix 
we prove a discrete-time version of a continuous-time result due to Datko which is required in 
section 4. 

2. The Discrete-Time LOR Problem-Feedback Solution and Approximation Theory 

In this section we briefly outline and summarize the infinite dimensional discrete-time 
linear-quadratic theory developed in [17] and the approximation results from [9]. Let X be a 

Hilbert space with inner product (v)x and corresponding induced norm M x . We consider the 
optimal control problem given by 

(<P) Find u = {u k } k _ 0 £ R m ) which minimizes the quadratic performance index 
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( 2 . 1 ) 


j ( u » x 0 ) = 2 (Qx k’ x k) x + u i Ru k 

k=0 

subject to the linear discrete-time control system 

(2.2) x k+1 = Tx k + Bu k , k = 0,1,2,... 

(2.3) x 0 e X 

where T, Q e 25 (X) with Q nonnegative self-adjoint, B e 25 (R m , X), and R is an m x m positive 
definite symmetric matrix. 

An input sequence u = {uj^} k _ 0 e °°> Rm ) * s called admissible for the initial data x Q if 

J(u;xq) < oo . An operator II e 25 (X) is called a solution to the algebraic Riccati equation 
corresponding to the plant defined by T,B,Q, and R if it satisfies 

(2.4) n =T*(II -ITB(R+ B*IlB)- 1 B*n)T + Q. 

We have the following theorem concerning the solution of problem ( V ). 

Theorem 2.1 There exists a nonnegative self-adjoint solution to the algebraic Riccati equation (2.4) 
if and only if there exists an admissible control for each x 0 8 X. If there exists an admissible 

control for each Xq e X, then the unique solution to problem (P) is given in linear state feedback 
form by 

u k = -Fx k , k = 0,1,2,... 

oo 

with the corresponding optimal trajectory x = {x k } k _Q given by 

x k+1 = Sx k , k = 0,1,2,... 

x 0 =x 0 

where F e 25 (X, R m ) and S e 25 (X) are given by 

(2.5) F = (R + B^B^B*!! T, 
and 

S = T - B F, 

andlle 25 (X) is the minimal nonnegative self-adjoint solution to (2.4). We have min J(u;x 0 ) = 
J(u;xq) = (IIx 0 ,Xq) x . If, in addition, any admissible control drives the state x k to zero 

asymptotically as k — ><=o (i.e. jim lx fc l x = 0; this would be true for example, if Q > 0) 

then the Riccati equation (2.4) admits a unique nonnegative self-adjoint solution II. If there 
exists an admissible control for each x Q e X, and Q > 0, then the spectral radius of S is less 
than one and S is uniformly exponentially stable. In particular if Q > 5 > 0 then I "S^l < Mr k , 
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k = 0,1,2,... with M = I n I / 8 and r = 1 - 6/ 1 n I < 1. 


In our discussion of the approximation theory below, we shall assume that the following 
hypothesis is satisfied. 

(HI) There exists an admissible control for each x 0 £ X, and any admissible control drives the 
state to zero asymptotically. 

We note that since F given by (2.5) is an element in Z (X, R m ) we have 
U k = -Fx k =-(f,x k ) x , k = 0,l,2,... 

where f = (fj, f 2 ,..., f m )^ with fj e X, j = 1,2 m, is referred to as the optimal functional 

feedback control gain. 

The approximation theory in [9] is based upon the finite dimensional approximation of the 
space X and the operators T,B, and Q. For each N = 1,2,... let X N be a finite dimensional 
subspace of X and let P N : X — > X N be the corresponding orthogonal projection of X onto X N . 
We take T N , Q N e £(X N ), e £(R m , X) and consider the following sequence of finite 
dimensional linear-quadratic regulator problems. 


(<P N ) Find u N = { u k } k=0 e i 2 (0» <*=; R m ) which minimizes the quadratic performance index 


yN/ N\ "V /n NN n. T _, 

1 (u;x 0 ) = 2u (Q + u k Ru k 


k=0 


subject to the linear discrete-time control system 


( 2 . 6 ) 


x k+1 = T^x* + b\, k = 0,1,2,. 


(2.7) 


N ..N 

X Q = P X Q £ X 


We require that the following hypothesis holds. 


N M 

(H2) For each N = 1,2,..., there exists an admissible control for every x 0 eX N and any 


N 


admissible control drives the state x k to zero, asymptotically as k — > °° . 


It then follows from Theorem 2.1 that for each N = 1,2,... there exists a unique, nonnegative 
self-adjoint solution fI N £ Z>(X N ) to the algebraic Riccati equation 


n N = (t n )* (n N - n N B N (R + (B^^nNB^-^BN)* n N )T N + q n . 
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The unique solution to problem (<P N ) is given by 


u? = k = 0,1,2,... 


with 


*N _ -o N x N 
X k+1 b X k ’ 


k = 0,1,2,...., Xq=P N x 0 , 


where 

fn = (R + (b n )* ri N B N )- 1 (B N )*n N T N 


and S = - B N F . We have min J N (u;Xq) = J N (u N ;x^) = (FI x^ ,x^ ) x , S has spectral 

radius less than one and is uniformly exponentially stable. Since F N e 35 (R m , X N ) it follows that 

u" = - (f 1 *,x£) x , k = 0,1,2,... 

where f = ,..., f m ) with f. eX , j = 1,2,..., m. 

The convergence theorem requires that the following hypothesis be satisfied. 

(H3) For each <J) e X, P N 4> -> <>, T N P N <> -> T(J», (T N )*P N <t> T*<J», and Q N P N (J) -> Q(J) as 
N -4 =«, and B N — > B in 35(R m , X) as N — > 

Theorem 2,2 Assume that hypotheses (HI) - (H3) hold. Suppose further that there exist positive 
constants Mj, M 2 , and r 2 , independent of N with r 2 < 1, for which 

(2.8) lIl N l < Mj 

and 

(2.9) l(S N ) k l < M 2 r 2 , k = 0,1,2 

Then for each $ e X, n N P N (J) — > n<}> and "S N P N (j) — > S<J), as N — » F N P N (j) — > F in 35 (X, R m ), 

— N — _H - -N - 

fj — > fj in X, j = l,2,...,m, and u k — > u k in R m and x k — > x k inX, for each 

k = 0,1,2,..., as N — » °o. 

3. Abstract Parabolic Systems with Unbounded Input 

Let H be a Hilbert space with inner product (•,•) and corresponding induced norm 1*1. Let V be 
another Hilbert space with inner product <-,-> and corresponding norm INI and assume that V is 
densely and continuously embedded in H with Icpl < pi||cp|| for cp e V. If we identify H with its 
dual, II*, it then follows that V <z H = H* c V* with H densely and continuously embedded in 
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V , the dual of V, endowed with the usual operator norm II- IK We have llcpll* < |i]cp| for (peH 
and Ikpll < |l 2 ||(p|| for (peV. We assume further that the embedding V c H is compact. 

Let a(',-):VxV-)C be a bounded, coercive, sesquilinear form on V. That is, there exist 
real constants a, y>0 and (3 for which 

Re a(cp,cp) > a||(p]| 2 - (3|cp| 2 , <p e V 
la(qw)| ^ Tllcpll INI, cp, V e V. 

The form a(y) defines an operator A £ £(V,V ) via 
(^9)(V) = C®P>V) = — a(cp,\|/), (p,\j/ e V 

where (•,•) in the above definition denotes the natural extension of the H inner product to the 
duality pairing between V and V*. If we define Dom(^) = [cp e V : jftp e H}, then Dom(^) = H 
and H : Dom (j?)cH-> H is the infinitesimal generator of an analytic semigroup { cT(t) : t > 0} 
of bounded linear operators on H with I <T (t)l < e®' 0 ^ 1 ,t > 0 (see [16]). In addition, it is shown 
in [16] that { <T(t) : t > 0} can be extended to an analytic semigroup on V , and in [2] that it can be 
restricted to an analytic semigroup on V. 

It is not difficult to argue (see [16]) that the H adjoint of Si, Si : Dom (Si ) c H —> H, is given 
by Si cp = ip, where is that element in H for which -a(0,cp) = (\j/,0), for all 0 e V, and 
Dom (Si ) consists of all those elements cp e V for which such a \j/ e H exists. The operator Si 
extends to an operator in £(V,V ) via 

(^*<p) (\p) = (A*<p, \j/) = -a(\|/,<p) s -a*(cp,\| /) , 

for cp, e V. It immediately follows that the sesquilinear form a*(y) :VxV-)C defined 
above satisfies 

Re a*(cp,cp) > a||cp|| 2 - p|cp| 2 , cp e V 
la*(<p»V)l ^ 7ll9ll IMI. <P. ¥ e V, 

and consequendy that Si is the infinitesimal generator of an analytic semigroup { cT (t) : t > 0} on 
V,H, andV* with l£f*(t)l < ,t > 0. It also follows that tT*(t) = £T (t)*, the H adjoint of 

«T (t), for all t > 0. 

We consider the continuous time control system 

(3.1) x(t) = j^x(t) + Bu(t), t>0 

(3.2) x(0) = x 0 

where £(V,V*) is as it was defined above, U3 e S(R m ,V*), u e L 2 (0, R m ), and 
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xq e H. We note that the fact that 73 is assumed to have range in V rather than H indicates that 
our framework will be able to handle certain classes of unbounded input - for example, certain 
types of boundary control. This will become clearer when we discuss examples below. We shall 
be concerned with the so called mild solution to the initial value problem (3.1), (3.2). The mild 
solution to the system (3.1), (3.2) is the function x e 1^(0,^ ;V) fl C(0, t f ; H) Pi H^O.tp V*) for 
any t f > 0 given by 


(3.3) x(t) =cT(t-s)x(s) + j* <T (t - o)73u(o)da, 0<s<t<T, 

S 

(3.4) x(0) = x Q , 

• • • • • jjc 

where the integral in (3.3) is interpreted as an integral in V . 

To derive the discrete-time system of the form (2.2), (2.3) corresponding to the system (3.1), 

(3.2) , we let x > 0 denote the length of the sampling interval, and consider piecewise constant 
(zero-order hold) controls of the form 

u(t) = uj, , t e [kx, (k+l)x), k = 0,1,2,.. 

where for each k, u^ is a constant vector in R m . Defining xj. = x(kx), k = 0,1,2,... , from 

(3.3) , (3.4) we obtain 

xj^ = Tx^ + Bu^, k = 0,1,2... 

Xq E H 

X 

where T e 35(H) and B e 35 (R m , H) are given by T = <T (x) and B = J 0" (t)73dt, respectively. 

o 

In setting up the LQR problem, in the performance index (2.1) we assumed that Q e 35(H) is 
nonnegative self-adjoint and that R is an mxm positive definite symmetric matrix. We make the 
following standing assumptions. 

(A) The pair {T,B} is uniformly exponentially stabilizable. That is, there exists an operator 

F e 35(H,R m ) for which the operator S = T - BF is uniformly exponentially stable; i.e. there 
exist positive constants M and r with r < 1 for which IS^I < Mr k , k = 0,1,2,... 

(B) There exists a 8 > 0 for which Q > 8 (i.e. (Qcp.cp) > 8[cpl 2 , cp E H). 

Lemma 3.1 If assumptions (A) and (B) hold then hypothesis (HI) is satisfied. 

Proof Let x q eH be given and set u^ = -Fx k , k = 0,1,2... where F is the operator in 
35(H,R m ) guaranteed to exist by assumption (A). Then x k = S k x 0 , u k = -FS k x 0 , k = 0,1,2... 
where S = T - BF E 35(H), and 

oo oo 

J(u; x 0 ) = (Qx k , x k ) + u^Ru k = ^ (QS k x Q , S k x Q ) + (FS k x 0 ) T RFS k x 0 

k=0 k=0 
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< J IQI IS k x/ + IRI IFI 2 IS k x/ < {|Q| + IRI IFI 2 } M 2 

k=0 


oo 

s 

k=0 


2k 


r 


{IQI + IRI IFI 2 }M 2 /l-r 2 <oo. 


Also, Q > 5 > 0 implies that any admissible control must drive the state x k to zero asymptotically 
as k — » oo , and the lemma is proved. 

4. Galerkin Approximation and Convergence 

For each N = 1,2,... let H N be a finite dimensional subspace of H with H N c V for all N. Let 

P N : H — > H n denote the corresponding orthogonal projection of H onto H N . We make the 
following assumption concerning the approximation properties of the subspaces H N . 

(C) For each cj>eV, lim IIP N <J> - <})ll = 0. 

Note that assumption (C) implies that P N cp — > cp in H as N — > °° for each cp e H and that 
||P N (p - (p|| < v[|(p|| for all (p e V and some v > 0 which does not depend upon cp or N. 

We use a standard Galerkin approach to define the operators e £(H N ). For cp N e H N , 

Let j^. N (p N = \|/ N where \p N is the unique element in H N guaranteed to exist by the Riesz 

Representation Theorem which satisfies -a(cp N ,0 N ) = (\)/ N ,0 N ) for all 0 N 8 H N . If we set 
tf N (t) = exp(j3. N t), t > 0, then { cf N (t) : t > 0} is a semigroup of bounded linear operators on 
H N with l<T N (t)l < e ^ , t > 0, N = 1,2,... . The adjoint of ^ N , (^ N )* 6 Z>(H N ), is given 
by (.T N )*cp N = \p N where \(/ N £ H N satisfies -a*(cp N ,0 N ) = (\|/ N ,0 N ) , 0 N bH n . It follows that 
<T N (t)* = exp((.a N )*t), t > 0. 


Lemma 4. 1 If assumption (C) holds, then for each cp e H we have cf N (t) P N cp — > <T (t)cp, and 
^r N (t)*P N cp -» r ( t )V as N — » oo, uniformly in t for t in bounded subintervals of [0,°°). 

Proof The result will follow from the Trotter-Kato semigroup approximation theorem (see [12]) 

once we have shown that (^ N - X) _1 P N cp -» (A- X ) _1 cp and ((;? N )* - A,) _1 P N cp -> (A* - X ) -1 cp 

as N — > oo for each cp e H and some X > (3. We argue that the first convergence holds only; the 
proof that the second holds as well is completely analogous. 

Let cp e H and set \j/ N = - X) _1 P N cp and \j/ = (A- X ) _1 cp for some fixed X > [3. (Note that X 

> P implies X e p (.3) fl p (.ft N ) for all N = 1,2,...). Then 
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ally* - P N y || 2 < Re a(V* - P N y, V* - P N y) + PI V* - P N \|/| 2 
= -Re((^ - X)(y N - P N \|/), y* - P N \|/) + (p - X) I y* - P N y | 2 

< -Re(P N cp, y* - P N \|/) + Re ((^ - X)(P N \/ - \|/), y N - P N \|/) + R e((JZ - I)\(/,y* - P N \|f) 

= Re(cp - P N (p,y* - P N y) - Re a(P N y - y,\jf N - P N y) - X Re(P N y - y/y 1 * - P N \j/) 

< i(p - p N cp|| y* - P N \|/ 1 + y||P N y - M\ II V* - p n v!I + MP*V - M 1^ - P N ¥l 

< KD-PV + e|y N -P N y| 2 + IIP N y-y]| 2 

+ y e ||y* - P N y || 2 + ^ IP N \i/-\j/f+ e|\/*-P N \(^ 2 

for any e > 0. Recalling that !0| < p||0|| for 0 e V, and choosing e > 0 so that a - e(2|i 2 + y) 
a/4, we obtain the estimate 

liy* - P N y [| 2 < (oce) _1 {|(p - P N cpl 2 + yll P N y - yll 2 + WP N y - yl 2 }, 

the right hand side of which tends to zero as N -» °° by assumption (C). It follows that 

ly 1 * - yl < ly 1 * - P N y| + IP N y - y| < P-Hy 1 * - P N y[| + IP N y - y] 

which tends to zero as N — » by the previous estimate and assumption (C), and the lemma is 
proved. 

Setting T n = $T n (t) 8 33 (H N ), Lemma 4.1 implies T N P N cp -> Tcp and (T N )*P N cp -> T*cp 
as N-^« for each 9 e H. 

Define ift N e 33 (R m , H N ) by <ft N v = y 1 * , where for v e R m , y 1 * is that element in 

which satisfies (again by the Riesz Representation Theorem) (Bv)(0 N ) = (y 1 *, 0 N ), for all 

0 N eH N . Define B N e33(R m ,H N ) by B N = J $T N (t)B N dt. 

0 
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Lemma 4,2 If assumption (C) holds, then B N — > B in £(R m , H) as N —> «> . 


Proof We assume (3 > 0. If p < 0, then P can be taken equal to zero in the argument that 

t t 


1 


follows. For t> 0, v 8 R'“, and N = 1,2,... set z(t) = | <r(s)tftvds andz N (t)= | «T* (s)(JV'v ds. 
Then z, z N e L 2 (0 ,t; V) ft C(0, x; H) fl H^O.x; V*), 




,N, N 


z(t) = MO +Bv, t > 0, z(0) = 0 
z N (t) = ^z N (t) + B N v , t>0, z N (0) = 0, 


and Bv = z(x) and B N v = z N (x). It follows that 

T 

a Je 4pt ilz N (t)-z(t)ll 2 dt 
0 

T 

< J e' 4(5t Re a(z N (t) - z(t), z N (t) - z(t)) + e 4pt p I z N (t) - z(t)l 2 dt 
0 

T 

= J-e^ReC^t) - MO, z N (t) - P N z(t)) + pe~ 4pt !z N (t) - z(t)l 2 dt 
o- 

T 

- J - e’ 4Pt Re(j ; lz N (t) - J3z(t), P N z(t) - z(t))dt 
o 


= J-e‘ 4Pt Re(z N (t)-z(t),z N (t) - P N z(t)) + p e' 4Pt lz N (t) - z(t)l 2 dt 
o 


+ J e’ 4Pt Re(.3z(t) - *z N (t), P N z(t) - z(t))dt 
o 
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+ J e 4(k Re(j3z(t) - ilz N (t), P N z(t) - z(t))dt 
o 

where in the last equality above we have used the fact that the definition of generalized derivative 
implies that 

X 

J (^{P N z(t)}-z(t),0 N (t))dt = O 
o 

for all 0 N e 1^(0, x; V). Continuing* we find that 

T 

a J e' 4pt llz N (t)-z(t)ll 2 dt 
o 

< J -e" 4pt Re(z N (t) - ^{p N z(t)},z N (t)-P N z(t))+ 2pe 4p ‘lz N (t) - P N z(t)l 2 dt 
o 

+ 2P J e' 4pt IP N z(t) - z(t)l 2 dt + J e^ pt Re(^z(t) - ^z N (t), P N z(t) - z(t))dt 
0 0 

= J‘ 2 ^ e' 4pt lz N (t) - P N z(t)l 2 dt + 2p JV 4pt iP N z(t) - z(t)l 2 dt 
o o 

T 

+ J e' 4Pt Re (^z(t) - ^z N (t), P N z(t) - z(t))dt 
0 

X 

< - j- e' 4pT lz N (x) - P N z(t)I 2 + 2pJ IP N z(t) - z(t)l 2 dt 

o 

x 7 x 

+ a J llz(t) - z N (t)ll 2 dt + j HP N z(t) - z(t)ll 2 dt, 

0 0 

and hence that 
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I *2 40x I 

Iz^x) - P N z(t)I 2 < 4p e 4Px J IP N z(t) - z(t)l 2 dt + - J llP N z(t) - z(t)ll 2 dt . 

o o 

Assumption (C), the remark immediately following it, the fact that z e 1^(0, ■ t;V) fl C (0,x;H), 

and the dominated convergence theorem therefore imply that IB N v - P N Bvl — > 0 as N — » °° for 

each v e R m . The triangle inequality, assumption (C) and the finite dimensionality of R m then 
immediately yield the desired result. 

If we define Q 1 ^ £ £(H N ) by Q N = P N Q then Q N is nonnegative (in fact positive, by 
assumption (B)), self-adjoint, and by assumption (C) satisfies Q N P N cp — > Qcp as N — > °° for each 
cp E H. This together with Lemmas 4.1 and 4.2 yield the following result. 

Lemma 4.3 For H N , P N , T N , B N , and Q N as defined above, assumption (C) implies that 
hypothesis (H3) is satisfied. 

To verify that hypothesis (H2) and the conditions of the convergence theorem, Theorem 2.2, are 

satisfied we shall require the fact that the operators B N : R m — » H N are uniformly bounded in 

3S(R m , V). Toward this end, we note that it can be argued (see [2]) that !l tT N (t)ll = !! <T N (t)II < CeP 1 , 
t > 0 where C is a positive constant which does not depend upon N. 

Lemma 4.4 If assumption (C) holds then there exists a constant Lj > 0 which does not depend on 
N for which IIB N vII < L^lvl, veR m ; that is the operators B N are uniformly bounded in 
£(R m ,V). 

Proof For veR m we have 

(XllB N v II 2 < Re a(B N v, B N v) + plB N vl 2 

= -Re(^ J $T N (t) (ft N vdt, B N v) + p(B N v, B N v) 
o 

= Re (^ N v - tr N (x) ?3 N v, B N v) + p(B N v, B N v) 

ll 



= Re ((ftv, B N v) - Re((Bv, B N v) + P(B N v, B N v) 

< llftvll* IIB N vll + Ce pT liBvIl* HB N vll + lp| ||B N vll* IIB N vll 

< lift *11 { 1 + Ce pT } Ivl IIB N vll.+ ipi|ilB N vl IIB N vll. 

Thus 

IIB N vII < (inB*ll/oc){l+ Ce Pt }lvl + (lp!p. / oc) IB N vl, 
which together with Lemma 4.3 yields the desired result. 

Lemma 4.5 Suppose that assumptions (A) and (C) hold and let F e £>(H, R m ) be the operator in 
assumption (A) which uniformly exponentially stabilizes the pair {T,B}. Then for all N 
sufficiently large, F uniformly exponentially stabilizes the pairs {T N ,B N }. That is, the pairs 
(T n ,B n } are uniformly exponentially stabilizable, uniformly in N for all N sufficiently large. 

Proof Let S^ = T N - B N F and let S = T - B N F. Then Lemma 4.3 implies that S N P N cp -> Sep and 

(S N )*P N cp — > S*cp asN— for each cpeH. It follows that there exists a constant Kj > 0 for 
which IS N I = l(S N )*l < K 1 and ISI = IS *1 < Kj , and consequently that the spectra of the operators 
S, S*, S N , and (S N ) are contained in the closed disc { z : Izl < Kq } in the complex plane. Let 

r\ > Kq. Then rj e p(S) 0 p(S N ) and since S N P N converges strongly to S as N — > oo, and the 

embedding V c H is compact, it follows (see [1]) that there exists a constant K 2 > 0 which does 

not depend on N for which I(rj - S N ) -1 P N i < K 2 . Therefore for cp £ H we have 

l(‘n-S N )- 1 P N q>-(Ti-S)- 1 <pl 

< l(n - S N ) _1 P N cp - P N (ri - s^cpl + l(P N - I)(tj - S) _1 cpl 

= i(ri - S N )V(S V - s)Cn - S)' 1 (pl + |(P N - m - s) Vl 

< KJ(S N P N - S)Cn - S) Vl + l(P N - I)(rj - S)'Vl 

z f 
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which by Lemma 4.3 tends to zero as N -» Analogously it can be shown that 

(H - -> (T) - S*)' 1 ^ as N -» °° for each tp 8 H and rj > Kj . 

We claim next that for some positive integer Nq we have 


(4.1) 


0 cr(S N 


N=N„ 


0(S ) c {z : Izl < e} 


for some e < 1. Suppose not. Then there exists a sequence of positive integers {Nj} J=1 


N: 


N: 


with Nj oo as j -> oa and sequences of complex numbers {A. J } J -_ 1 and elements in H J , 

N- ~ N- N- N- 

{cp J } j =1 for which [ X J , cp J ] is an eigenvalue / eigenvector pair for the operator S J with 

N- N- N- N- ~ 

IX, J I > 1 - 1/j andltp J 1 = 1. Now IX. J l<Kj. Therefore {X, J }j=i must admit a convergent 

N. 

subsequence {X, }“ =1 with limit i satisfying 1^1 >1. For convenience we re-index and 

say lim X N = 1. We claim that 1 is an eigenvalue of S. Indeed, for rj > K. we have 

(t) - X N )tp N = (r| - S N )cp N , or 


(4.2) (p N = (t) - S N r ! (rj - X N ) 9 N . 

The estimate (see [16]) ll<T N (t)\(/ N ll < — hj/ N l, t>0, y N e H N , for some positive constant 

■ft 

L 2 which does not depend on N, and Lemma 4.4 imply 


IX, N | llcp N ll = llX, IN cp N ll = IIS N (p N ll < IISr IN ('t)(p N ll + IIB N F(p ] 


i n „N m _ llc N n, 


N ^,.N,i , i,r,N„ Nil 


< hi |(pN| + T IFI |(pN| = hi + L.IFI 

Jx v 1 v Jx 1 


Recalling that X N X and | X | > 1, it follows that the sequence {cp N } lies in a bounded subset 


of V. The assumption that the embedding VcH is compact then implies that {cp N } admits an 


N 


H-convergent subsequence {(p J } with limit tpeH. Once again we re-index and assume that 

lim l(p N - cb| = 0. Then for any \j/eH, (4.2) yields 
N— 

(On - 01 - (S N )*) _1 P N X 1 /) = (cp N , p\) = (cp N ,v). 
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Taking the limit as N->»we obtain 


(Cn - £)$, Oi - s*)'V) = ($, ¥)> 

or Sep = £ cp; i.e. [ X, cp] is an eigenvalue/eigenvector pair for the operator S. But this is a 

contradiction since l£| > 1 andS=T-BF is assumed to be uniformly exponentially stable. It 
follows that (4.1) must hold for some positive integer Nq sufficiently large. 

Now for N > Nq, let x N e H N and set x k = (S N )^x N . If we denote the z-transform of the 
sequence {x k } k _ 0 by x N (z), then x N (z) = (I - z" 1 S N )- 1 x N , and 

x t = 2^ J 

where the integral in the above expression is around a closed contour in the complex plane which 
contains the spectrum of the operator S N in its interior. If we choose the contour to be the circle 

z = re i0 , 0 < 0 < 2k with e < r < 1, we obtain lx k [ < K 2 r k+1 lx N l, k = 0,1,2,..., or 

|(S N ) k I * Mr* k = 0,1,2,... 
with M = K 2 r and the proof is complete. 

Theorem 4.1 If assumptions (A), (B), and (C) hold, then the conditions of Theorem 2.2 are 
satisfied and therefore the convergence results stated as the conclusions of that theorem are valid for 
the class of problems considered in section 3. 

Proof . Lemma 3.1 implies that hypothesis (HI) is satisfied and Lemma 4.3 implies that hypothesis 

(H3) holds. Assumption (B) and the definition of Q N yields (Q N cp N ,cp N ) = (P N Q(p N , tp N ) = 

(Qcp N , cp N ) > S[(p N l“, cp N e H n . This together with Lemma 4.5 and the same arguments used to 
prove Lemma 3.1 imply that hypothesis (H2) is satisfied. Thus we need only to verify (2.8) and 
(2.9). Toward this end we note that 

(4.3) in N | = sup {(n N cp N , cp N ) : cp N e H N , lcp N l = 1} = sup J N (0 N ;cp N ) 

191=1 

where u N is the optimal control corresponding to the initial data cp N . For cp N e H N with 
|cp N l=l define u N ={u k ) k=0 byu k =-Fx k , k = 0,1,2,... where F e £(H, R m ) is the 

N M 

operator from assumption (A) and {x k } k=0 is given by (2.6), (2.7) with x Q = cp . Then for N > 
N 0 with S n = T-B n F, Lemma 4.5 implies 
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)VV)S Au N ;<P N ) = £ (Q N (S N M (S N )V) + (F(S N ) k <p N ) T RF(S N )V 

k=0 

£ £|Q N II(S N )Vi 2 + IRI IFI 2 l(S N ) k <p N | 2 < (IQ| + IRI IFI 2 )M 2 kp N l 2 £ 

= {IQI + IRMF| 2 }M 2 /(l-r 2 ) =M l 
which together with (4.3) establishes (2.8). 

Now in N l uniformly bounded implies that the operators = T N - are uniformly 

bounded. Thus 

(4.4) l(S N ) k l <s k , k = 0,1,2... 

for some positive constant s which does not depend on N. Also, for tp N e 

(4.5) £ (Q N (S N ) k cp N , (S N ) k tp N ) < (nV>9 N ) ^ M^l 2 . 

k=0 

Since > 5>0, (4.4) and (4.5) together with a discrete-time version of a result due to 
Datko [5] (see the Appendix) establishes (2.9) for some M 2 ,r 2 > 0 with r 2 < 1 and the theorem is 
proved. 

Finally we note that once a basis has been selected for the finite dimensional subspace H N , the 
matrix representations for the operators T N ,B N and Q N , which are required to solve the 

N xr 

approximating optimal control problem (IP ), are easily computed. Indeed, if {<Pj }j =1 

denotes a basis for H N , then [^L N ] = - (cpj , cpj ) _1 a(cp^ , cp^ ) T with [T^] = exp (t[^ N ]). 
Similarly we find [ft N ] = (cpj , cpj ) -1 ( (fce k , cp^ ) T where denotes the k* standard unit 

X 

vector in R m , [B N ] = | exp (t[#])[B N ]dt, and [Q N ] = (cp*, (p^) -1 (Qcp™, cp^). A more 

o 

complete and detailed treatment of the computational and implementational aspects of solving the 
finite dimensional approximating linear-quadratic regulator problems (tP N ) can be found in [9]. 

5. Examples 


We describe a generic class of boundary control systems and some specific examples to which 
the abstract approximation framework and theory we have developed above applies. Let H be a 
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Hilbert space with inner product (•,•) and corresponding induced norm 1-1, and let W be another 

Hilbert space with inner product and norm [•,•] and lll-lll respectively. Let <•,•> be another 

inner product defined on WxW and let the corresponding induced norm he denoted by INI. We 

assume that WcH and that there exist positive constants |i and p for which 

Icpl <|ill(pll < plllcplll, forallcpeW. 

We consider the generic boundary control system given by 


(5.1) 

x(t) = Ax(t), 

t>0 

(5.2) 

Tx(t) = u(t). 

t>0 

(5.3) 

x(0) = Xq 



where x(t)eH, u(t)eR m , x 0 eH, Ae£(W,H) and T 8 £(W, R m ) is assumed to be 

surjective. Let Wq denote the null space of T, 71 (T). That is, Wq = {cpeW : Ftp = 0 } and define 

Y to be the completion of Wq with respect to the norm INI. It follows that V is a Hilbert space 

with inner product and corresponding induced norm ll ||. We shall require the following 
assumptions: 

(1) V is dense in H 

(2) The embedding VcH is compact 

(3) W c V. 

Then choosing H as our pivot space, we have V c H c V with the embeddings dense and 
continuous. 

Define the bounded sesquilinear form o(-,-) : W x H — » C by a((p,\|/) = -(Arp/vp) for (p e W 
and \|/ e H, and assume that there exist real constants a,P,y with a,y > 0, for which 

(4) Reo((p,<p) > otllcpll 2 - Plcpj 2 , (p e W 0 

(5) lo((p,\j/)| < y||(p|| ||\j/||, cp,\j/eW 0 . 

Define the sesquilinear form a(-,-) : W 0 x W 0 —> C to be the restriction of the form a(v) to 
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W 0 x W 0 . It follows from density and continuity (i.e. assumption (5) above) that a(y) admits a 
unique extension to a bounded coercive form on V x V. Thus assumptions (4) and (5) continue to 
hold with a(y) replaced by a(y) and the space W Q replaced by V, and consequently we may 
define an operator J48 3o(V,V )via 

C®p)(\|0 =(Jk p,y) = -a(<p,\|/), cp,\|/ £ V. 

If we restrict !A. to Dom (A) = {cp e V : Jftp 8 H} then A : Dom( J$) a H — > H is densely defined 
and is the infinitesimal generator of an analytic semigroup { <T(t) : t > 0} of bounded linear 
operators on H with I <T (t)l < e(P -a/ tiH t > 0. Also, { O' (t) : t > 0} admits an extension and a 

• • ;fc 

restriction to an analytic semigroup of bounded linear operators on V and V, respectively. 

Recall that P was assumed to be suijective and let r + £ £ (R m , W) denote a fixed but arbitrary 
right inverse of T. That is rr + u = u for u e R m . Define B e 3S(V,V*) by 
(5.4) (fau)(<p) = (Bu,cp) = air 1 " u,(p) - c(r + u,cp), 

for u 8 R m and cp e V. We note that the operator "53 is indeed well defined. For if T 1 and T 2 

are two distinct right inverses of T and we let and denote the corresponding operators 

+ + + + 

defined as in (5.4), then for any u e R m and (peV we have TjU - r 2 u = (r j - r 2 )u e W 0 and 
(ITi x u)((p) -0ft 2 u)(cp) = a((ri - r 2 )u,(p) -ct((Ti - T 2 )u, cp) 

= o((Ti - r 2 )u, cp)-c((ri - r 2 )u,cp) = o, 

or ^3 = B 2 . 

Following Curtain and Salmon [4], for any t f > 0 we call x e L 2 (0, t f ; V) fl C(0, t f ; H ) fl 
H^O.tj; V*) the weak solution to the boundary control system (5.1) - (5.3) if it is the unique mild 
solution to the initial value problem 

x(t) =31x(t) + fl3u(t), t>0 
x(0) = x 0 . 

That is, if 

t 

x(t) = tT (t - s)x(s) + j* <T(t - a)1!3u(a)dcj, 0 < s < t < t f . 
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x(s) = X 0 . 

For the discrete-time system, with x denoting the length of the sampling interval, we have 

X 

T = cT(x) e £(H) and B = j*<T (t)<Rdt e £(R m , H). From (5.4) together with the definitions 

o 

of the forms g(v) and a(v), we obtain 

XX XX 

Bu = jV(t)ftudt = jV(t)(A - J3)r + udt = J^T(t)Ar + udt- J^tT(t)r + udt 
oo oo 

X 

= (I - tT (x))r + u + J ?T (t)Ar + udt . 
o 

Similarly, for the approximating input operators B N we find 

X 

N N N + f N N + 

B u = (I-tr (x))F T u + KT (t)P AT udt. 

o 

If F*" is chosen so that <Tl(F i ') c 71>(A), the discrete-time input operators B and B N take on the 

particularly simple forms B = a - = (I - T)r + and B N = (I -$T n (t))P n r+= (I - T N ) P N I^. 

We discuss two specific examples of boundary control systems of the general form that we 
have just described, and outline the arguments necessary to verify that assumptions (A) - (C) are 
satisfied. 

Example 5.1 We consider a system whose dynamics are described by the one dimensional heat 
equation 

(5.5) (t/q) = a Q (t,r|), t > 0, 0<T| < 1 

where ag > 0, with a homogeneous Dirichlet boundary condition at T| = 0, 

(5.6) w(t,0) = 0, t>0 

and Neumann boundary control at iq = 1, 

(5.7) a Q ^ (t,l) = u(t), t>0. 

The initial conditions are assumed to be of the form 

(5.8) w(0,ri) = w 0 (r|) , 0 < ri < 1 
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where w Q e 1^(0, 1) is given. 

To put the system (5.5) - (5.8) in the form of (5.1) - (5.3) we take H = I^O.l), and 
W = H 2 (0,l)nHi (0,1) where (0,1) = {cpeHHO.l) : cp(0) =0}. The inner products 
(•,•)» <•»•>, and [•,•] are taken to be 


l 


(<P,V) = f w. 

(p,\)/eH 

V 

0 


i 

/• 


<cp,\|/> = DcpD\j/, 

9 ,\j/eW 

V 

0 


l l 

* m 

i 

* . _ 

[9,\|/] = 9V + 

D(pD\(/ + D cpD \\r, 

J 

J J 

0 0 

0 


The operators A E £(W,H) and T e ^(W.R 1 ) are given by Atp = agD 2 (p and Tcp = agD(p(l), for 
cpeW. We then have W 0 = %(T)= {(peH 2 (0,l): cp(0) = Dcp(l) = 0}, a( v ) : W x H C 
given by a(cp,\}/) = -(a 0 D 2 (p,\|/), cp 8 W, \\r e H, V = H[ (0,1), V* c H -1 (0,1), and a( v ) : 

V x V -» (D given by a(cp,\j/) = (aoDcp,D\j/), (p ,\j/ e V. The operator Az £ (V,V*) takes the form 
(^cp)(\(/) = - (aQD(p,D\j/) for cp,\j/ e V. It follows that A : Dom(^) cH-)H is given by j^cp = 
a 0 D 2 (p for (p e Dom(^) = Wq and is self-adjoint. We have 

a(cp,(p) >a 0 ll(p|| 2 , cpeV 
la((p,\|/) I £ a 0 Mcp[[ ||\|/||, 9>Y e V, 

and therefore that (^cp,cp) < - a 0 lcp| for 9 £ Dom(JI). Thus A is the infinitesimal generator of 
an analytic semigroup { <T(t) : t > 0} of bounded, self-adjoint linear operators on H (and V and 

5jC 

V ) which is uniformly exponentially stable; that is 

(5.9) IST(t)l < e‘ aQt , t>0. 

We take F 1- e S(R ! ,W) to be (I^u) (rj ) = (u/aQ)T] , 0 < tj <1 and note that with this choice 
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of r^we have BfU 4 ’) c ri(A). We note further that W is indeed a subset of V and thus 
B e ^(R^V*) as defined in (5.4) is given by (Bu)(cp) = u(p(l), for u e R 1 and cp e V. The 
discrete-time input operator B e S( R 1 ,''/*) will take the form Bu = (I - 0" (x)) T^u. We have 
x(t) = w(t,-), t > 0 and xq = wq . 

The space V is densely, continuously, and compactly embedded in H, it is clear from (5.9) that 
assumption (A) is satisfied and if, for example , we choose Qe 35 (H) to be Qcp = qcp where 

q £ £^,(0,1) is such that qCq) > 0, a.e. T| £ [0.1], then assumption (B) will be satisfied as well. 
With regard to approximation, a linear spline based scheme is one for which assumption (C) can be 

shown to hold. For each N= 1,2,... let {cp^} j_ 0 denote the usual N + 1 linear B-splines (i.e. 

"hat" functions) defined on [0,1] with respect to the uniform mesh [0, 1/N, 2/N,...,l }. Let H N be 

N xi N N N 

the N-dimensional subspace of V = (0,1) given by H = span {(pj } j=] , and let P : H — > H 

denote the corresponding orthogonal projection of 1^(0, 1) onto H N . Using the approximation 

properties of interpolatory splines (see [15]), the Schmidt inequality (see [16]), and the variational 

N 

properties of the orthogonal projection, it is not difficult to argue that lim IIP (p - cpll = 0 for each 

cp £ Hl( 0,1) and consequently that assumption (C) is satisfied. We note that other commonly used 
finite element methods for the heat equation can be shown to lead to approximation schemes for 
which assumption (C) holds, including for example, modal and spectral methods. 

For the system which has just been discussed, the optimal functional feedback control gain 

takes the form of a function f £ 1^(0, 1) with the optimal control given by 

l 

u k = - J f(t|)w(kx,ri)dri, k = 0,1, 

o 

The approximating optimal gain is a function fiN e h n with 

l 

u k =-Jf (fi)w N (kx,ri)dri, k = 0,1,2... 
o 

where w is the system state which results from the input u = { u k } k=0 • Our theory yields 
_N _ 

that lim If - f I = 0 and numerical studies that we have earned out and reported on elsewhere 
N— 

(see, for example, [7], [9], [10], [13]) substantiate our theoretical findings. 
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Example 5.2 In this example we consider the control of the transverse vibration of a cantilevered 
Euler-Bemoulli beam with Kelvin-Voigt viscoelastic damping and shear boundary control at the 
free end. The dynamics of such a system are described by 

(5.10) % (t,T|) + c 0 4*- (t,rj) + a 0 (t,ri) = 0, t > 0, 0 < T] < 1 

9t dry 01 dr) 

(5.11) w(t,0) = ^ (t,0) = 0, t>0 

or\ 

(5.12) c 0 i 3g-(t,l) + a Q pL (t,l) = 0, t > 0 

0 9rj 2 ot 0 3 tj 2 

(5.13) - c 0 (t,l) - a (t,l) = u(t), t > 0 

0 0ri 3 ot 0 dry 1 

(5.14) w(0/n) = w 0 (ri), 4^ (0,ti) = Wi(ti), 0 < ri < 1, 

where a Q , c Q > 0, and w Q e H 2 (0,1) = {cp e H 2 (0,1): (p(0) = D(p(0) = 0}, and Wj e L 2 (0,1) 
are given. 

To put the boundary control system (5.10) - (5.14) in the form of (5.1) - (5.3), we let 

2 

H = H l ( 0,1) x L 2 (0,1) with inner product 

t t 

((qytl/j), ((p 2 >V 2 )) = a o J D \ D \ + J ^1^2 > 

0 0 

and let W = { (cp,\|/) e H : c Q \j/ + a Q (p e H 4 (0,1), \j/(0) = D\|/(0) = 0, CqD^I) + a Q D 2 (p(l) = 0} 
endowed with the inner product [y] taken to be the standard inner product on H 4 (0,1). We 
define the inner product <• ,•> on W by 

l l 

<(<Pi>Yi)> (<P 2 >V 2 )> = a 0 J d 2 ( Pi d 2 ( P 2 + J D Vi D V 2 • 

o o 

i ^ ^ 

The operators A e 3S(W,H) and T e 35(W,R ) are given by A(cp,\|/) = (\|/, - c Q D \\r - a Q D (p) 
and r((p,\j/) = - c 0 D 3 \]/(l) - a Q D 3 (p(l), for (<p,\j/) e W. It follows that W Q = {(cp,\|/) e W: - 
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c 0 D 3 vj/(1) - a 0 D 3 cp(l) = 0} and that V = H^O.l) x H \ (0,1). The form o( v ) : W x H -» ID 
is given by 


l 


l 


l 


o((9»V). ( 0 >X)) = - a 0 1 D 2 \]/D 2 9 + a Q | D 4 cpx + c Q J D 4 \j/X 

0 0 0 

for ( cp,\j/) e W and (0,%) e H, and the form a(v) : V x V — » (E is given by 


l l l 

a((qVYi). (<P 2 .V 2 )) = - a 0 J D VjD 2 q> 2 + a 0 J D^D 2 ^ + c Q j* dVjDV 2 > 


for (<p.,\|/.) eV, i = 1,2. 

It is not difficult to show that 


a((qw)> (<P>V)) ^ c 0 ll((p,v)|| 2 - c 0 l((p,\|/)| 2 
la((cp,V), (9,X))I ^ 72max(l,c 0 ) !!(<p,\|/)|| ||(0,x)|| 

for ( cp,\j/) , (0,%) e V. With the operator j%£ £ (V, V*) defined by 


(*(qw)) ( 0 >X) = -a ((<p.V), ( 0 ,X)) 

for ( (p,\j/) , (0,x) e V, we have A : Dom(j?) c H -> H given by J3((p,\|/) = A(cp,\|/) for 
( (p,\{/) eDom(j?)=W 0 . The operator j?: Dom(j?)c= H —> H is densely defined and the 

infinitesimal generator of an analytic semigroup of contractions { <T(t) : t > 0} on H. Moreover, it 
can be shown (see [6]) that the semigroup { tT(t) : t > 0} is in fact uniformly exponentially stable. 
Let cp 0 be the cubic polynomial which satisfies the interpolator conditions 

(p 0 (0) = Dcp 0 (0) = D 2 ip 0 (l) = 0, D 3 9 q (1) = -l/a^. We choose T + e 35(R ! ,W) as T + u = (cp Q u,0), 

u e R 1 . Once again we note that fl(r + ) c r](A), and , since WcV, the operator B as defined 

by (5.4) is an element in £ (R^V*) with (Bu)(tp,\|/) = u\)/(l), ueR 1 , ( (p,\)/) eV. The 

discrete- time state transition operator is given by T = *T(t) e 35(H) and the discrete- time input 

operator is given by B = (I - T) r + e £ (R 1 , H) . 

The uniform exponential stability of the open-loop semigroup { cf(t) : t > 0} implies that 
assumption (A) is satisfied. In flexible structure control problems such as the one we have just 
described, the state penalization operator Q in the performance index is frequently chosen to be the 
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identity (i.e. the energy of the system is to be driven to zero). Thus when this is the case, we have 
that assumption (B) is also satisfied. 

N N+l 

For an approximation scheme, for each N = 1,2,... we let {B^ }^ = 1 denote the standard 

cubic B-spline functions defined on the interval [0,1] with respect to the uniform 

mesh {0, 1/N, 2/N,...,l}. Let denote the modified cubic B-splines which satisfy 

(3^(0) = Dp^(0) = 0, j = 1,2,...,N+1 (that is, B* = B* - 2B* - 2B N p (3* = B* j = 2,3,...,N+1), 
and set jf = ((3*0), j = 1,2,..., N+l, ft* +1+j = (0,(3*), j = 1,2,..., N+l. Let 

N 

N a 1 ' 9N+? N N N 

H = span {[3. } c V and let P : H — > H denote the orthogonal projection of H onto H . 

Once again elementary properties of spline functions (see [15]) and standard techniques from the 

N 

theory of approximation (see [14]) can be used to argue that lim^ ||P (cp,\|/) - ((p,\|/)|| = 0 for 

each (cp,y) e V, and consequently that assumption (C) is satisfied. 

The optimal functional feedback control gains are of the form f = (fj, f 2 ) £ H (i.e. 
fje H 2 (0,1), f^e L 2 (0,1)) with 

t l 

u k = -a Q | D 2 f 1 ('q)D 2 w(kx,Ti)dii - j* f 2 (ri)w(kT,Ti)dcp, k = 0,1,2,... 
o 0 

(note that in this example we have x(t) = (w(t,-), w(t,-)), t > 0). The approximating 

functional feedback control gains are of the form f = (fj , f^) E H* and our theory yields 
-N _ 2 -N - 

lim f, = f. in H and lim f„ = L in L„. Numerical studies for a flexible structure 
example such as the one we have described above were reported on in [8] and [13]. 

Finally we note that both of the examples that were treated above admit generalization to higher 
dimensions with the approximation subspaces H N being formed via tensor products of one 
dimensional elements. Also, it is worth noting that not every control system involving unbounded 
input that we might formulate will conform to our general framework. For example, if we replace 
(5.7) with Dirichlet boundary control, such a system results. In this case we have that W is not a 

subset of V and consequently the continuous-time input operator <J3 will have range in some space 
larger than V . More specifically ^3,(13) will be contained in the dual of Dom(^[ ) (see [4], [8], 

[10], [13]), and the range of the discrete-time input operator, (B>(<ft) will be a subset of W rather 
than V. While it is in fact possible to apply the theory outlined in section 2 to this particular 
example (see [9]), the question as to whether or not the results in section 3 and 4 can be extended 
so as to be able to handle this more general class of systems with unbounded input, at present, 
remains open. 
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6. Summary and Concluding Remarks 

We have developed a rather complete abstract approximation framework for discrete-time LQR 
problems for parabolic systems with unbounded input. Our results are in some sense analogous to 
those found in [3] for continuous time systems but also represent an extension in that the theory 
developed here can handle at least some, but certainly not all, forms of unbounded 
(continuous- time) input. Requiring only that (1) the infinite dimensional open-loop discrete-time 
system be uniformly exponentially stabilizable, (2) that the state penalization operator in the 
quadratic performance index be positive definite, and (3) that the usual Galerkin hypothesis on the 
finite dimensional approximation subspaces (i.e. the strong V-norm convergence of the 
corresponding orthogonal projections) be satisfied, we are able to apply the theory developed in [9] 
to obtain norm convergence of the approximating optimal feedback gains. 

In addition to the question raised at the end of the previous section, other related open problems 
remain. For example, in regard to the condition (1) above, one might ask if the uniform 
exponential stabilizability of the open-loop continuous-time system implies uniform exponential 
stabilizability of the corresponding discrete-time system. In a finite dimensional setting, this notion 
is treated in a paper by Hautus [1 1]. It is shown there that the answer to this question is yes for all 
but a finite number of sampling interval lengths z. The extension of this result to the most general 
of infinite dimensional systems is not immediately clear, and is worthy of some consideration. 

The extension of our theory to the LQG problem for infinite dimensional discrete-time systems 
with both unbounded input and output in the spirit of the treatment of this problem in [10] should 
also be looked at. In particular, since the semigroup { iT(t) : t £ 0} is analytic, the discrete- time 

input operator B is in fact an element in £(R m ,V). Therefore, if the output operator is bounded 
on V, one might attempt to develop a general theory wherein the problem is formulated in the state 
space V rather than H. However, without additional structure, for example, the operator A 
self-adjoint, such results may be difficult to obtain. 
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x: A Discrete-Time Version of a Result of Datko's 


Let X be a Banach space with norm M x and suppose Se3j(X). Suppose further that there 
exist positive constants Mj, M 2 and rj with r t > 1 for which IS k l < 1 , k = 0,1,2... and 


ca 

(A.l) ^ IS k xl x < M 2 lxl x , x e X. 


Then there exist positive constants M and r with r < 1 such that IS k l < Mr 151 , k = 0,1,2,... . 

Proof . Following Datko's proof in the continuous time case we first show that there exists an 
M 3 > 0 such that IS k xl x ^ M^lxl x , k = 0,1,2,..., xeX. Suppose not. Then there exist 
sequences {(k n , k n )} and {x n } with 0 < k n < 3c n < °° and lx n l = 1 such that 


(A.2) 


Mjn < IS^x L < M.r kn ' kn 

1 n A 11 


Then for k e [k n , k n ] an integer, we have 


which implies 


k n -k k-k n V .k k-k_ 

Mjn < IS S x n l x < M/f IS x n l x 


- (kn-k) 


nr, - IS x n l x . 


Squaring both sides of the last inequality and summing from k = k n to k = k n - 1, we obtain 


I -V, 


s X 


X IS K ' K ”x n l 2 s X S M 2 Ix "' 2 - ^ 
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Thus 


Vi 


X 


2_-2<V« = n 2 r ‘ V, 


VI 


X? 


n 2 (1 - 1 " 2 ^ n ~V^ 


k=k n 


k=k tl 






( l_ r - 2 <W) < 


..V^n 


Now (A.2) implies n < y n , or k n - k n > log^n. Therefore for some n Q sufficiently 

large we have 1 - r 2 ^ n ^ > y for all n > n Q . Thus for all n > n Q we have n 2 /2(r^ - 1) < M 2 
which is a contradicton. 

We next argue that given any 8 > 0 with 8 < 1 there exist a K e 2 with 0 < K< 1/e 2 
,K 


for which IS xl x < J M 2 e Ixl x for all x e X. The bound (A.l) implies that for each x * 0 
there exists aK = K(x,e) such that 

IS xl x < J IV^ 8 lxl x and £ lxl x < IS k xl x , 0 < k < K. 


This implies that 

K-l 

M^lxl 2 K < IS k xl x < M 2 !x£ 
k=o 

which in turn implies that e 2 K <1 or K = K(e) < 1/ e 2 . Thus K may be chosen independently of 
xeX and the claim is established. 

Now let 8 > 0 be given and choose e > 0 so small that ^M 2 e < 1 and M 3 ^/M 2 e < 5. 

Then for k > K(e) we have k = K , i > 0 and 

IS k xl x = IS K+ ^xl x = IS 4 S K xl x < MjIS K xl x < elxl x <8lxl x 

for any x e X. Therefore S is uniformly asymptotically stable. 

Since S is uniformly asymptotically stable there exists a K > 0, an integer, such that 
k 1 

IS I < y , k > K. Now any k = 0,1,2,... can be written as k = nK + m where n is a natural 
number and m is an integer such that 0 < m < K. Then for k = 0,1,2.. 

!S k | = |s nK+m | < IS m l IS^I < IS m IIS K | n <M 1 i^ 1 (l/2) n = M 1 e mlnri e’ nIn2 
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(mlnr 1+ (m/K)ln2) ( k/K )ln2 Klnr + In 2 _ (k/K )l n 2 

M x e e < M t e e 


k_ 

= 2M lt f (i) K = Ml k 


where M = 2M x r^ and r = 1/^2 < 1. 
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